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ABSTRACT 

In this work we propose a novel method for testing the vaUdity of the fiducial ACDM cos- 
mology by measuring the cumulative distribution function of the most massive haloes in a 
sample of subvolumes of identical size tiled on the sky at a fixed redshift. The fact that the 
most massive clusters probe the high-mass tail of the mass function, where the difference 
between ACDM and alternative cosmological models is strongest, makes our method particu- 
larly interesting as a cosmological probe. We utilise general extreme value statistics (GEV) to 
obtain a cumulative distribution function of the most massive objects in a given volume. We 
sample this distribution function according to the number of patches covered by the survey 
area for a range of different "test cosmologies" and for differently accurate mass estimations 
of the haloes. By fitting this sample with the GEV distribution function, we can study which 
parameters are the most sensitive with respect to the test cosmologies. We find that the peak of 
the probability distribution function of the most massive halo is well suited to test the validity 
of the fiducial ACDM model, once we are able to establish a sufficiently complete large-area 
survey with Mi™ ^ lO''*-^ Mq /i"' (Mii„ ^ 10^'^ M^h-^) at redshifts above z = 1 (z = 1.5). 
Being of cumulative nature the proposed measure is robust and an accuracy of 20 - 30% in 
the cluster masses would be sufficient to test for alternative models. Since one only needs the 
most massive system in each angular patch, this method would be ideally suited as a first fast 
consistency check before going into a more complex statistical analysis of the observed halo 
sample. 

Key words: Galaxies; clusters; general - Cosmology; miscellaneous - Methods; statistical 



1 INTRODUCTION 

Recently the study of the most massive galaxy clusters in the ob 
servable Universe saw an increased interest (IMantz et al.l 2008 



Cavon et al.| |201C|; IhoIz & Perlmutt"eil l2010l: iMantz et al. 



201C 



Baldi & Pettorindl201lL iHovle et alj201ll : iMortonson et alj201ll) . 
which was mainly initiated by the discovery of the very mas- 
sive high-redshift cluster XMMU J2235 .32557 at z = 1.4 with 
M200 = (7.3 ± 1.3) X lO'-* Mo ( MuUis et al.ll2005l ; iRosati et all 
I2OO9I ; IJee et ^|2009|) . Those studies mainly concentrated on the 
consistency of the presence of extremely massive clusters at in- 
termediate and high redshifts with the ACDM concordance model. 
Particular attention was also given to the impact o f non-Gaussianity 
on the high mass end of cosmological structures jjimenez & Verdel 



2009l;ICav6n et afbOldJEnqvist et alj201ll ; |ParaniaDe et al.l201ll : 
lOb 



Sartoris et al.ll201C 

Moreover, recently iDavis et al] (l201l[) applied genera l extreme 
value statistics (GEV) (see e.g. iFisher & TippetlllTgish : iGumbeil 
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( Il958h : IColesI ( l200lh ) to study the probability distribu t ion of the 
most massive halo in a given volume and Colombi et al.l ( I2OIII) ap- 
plied GEV to the statistics of Gaussian random fields. Apart from 
this, GEV, b eing relatively wide spread in the environmental (see, 
e.g.. lKatz & B rown ( 1992), Katz et al. (2002)) and the financial sci- 
ences (see, e.g., limbrechts & Schmidll ( il994l) ), has seen very few 
applic ations in the framework of astrophysics. iBhavsar & Barrowl 
( Il985l) . for i nstance, stud ied the statistics of the brightest cluster 
galaxies and IColesI ( Il988[) applied GEV on the temperature max- 
ima in the CMB. 

In this work, however, we are not interested in the single most 
massive halo in the Universe, but in recovering the cumulative dis- 
tribution function (CDF) of the most massive halo in subvolumes 
by fitting a CDF obtained from GEV and to study its possible dis- 
cerning power for testing different cosmological models. We will 
present an attractive, simple and robust method for model testing 
that c ould be applied to fu ture wide-fiel d cluster surveys lik e EU- 
CLID jLaureiis et al.l2009l) or eROSITA dCappelluti et alj2011i) for 
instance. 

The plan of this paper is as follows. In Sect. |2] we briefly intro- 
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duce the applicati on of GEV on massive clusters as discussed by 
lOavis et aUfeOllh . followed by an introduction to our idea of mea- 
suring the underlying CDF for massive clusters in Sect.|3] In Sect.|4] 
we study the parameters of the GEV distribution for several cos- 
mological test models. First, for the case of a small fixed comoving 
volume in Subsect. 14.11 and then we extend our method to subvol- 
umes of arbitrary depth in redshift space in Subsect. 14.21 In order 
to study the usability of our method, we utilise in Sect. |5] an in- 
verse sampling technique to create observed samples of the under- 
lying CDF distributions for the different cosmologies and recover 
them by fitting a GEV cumulative distribution function to them. 
Section|6]discusses the robustness of the method with respect to in- 
accuracies of the mass estimates. We summarise our findings in the 
conclusions in Sect.[7]and briefly review the aspects of GEV as it is 
necessary for this work including a list of the most useful relations 
in the AppendiceslAlandlB] 



2 THE GEV STATISTIC 

This section briefly introduces the relations of the GEV relevant 
for this work. A mor e detailed di s cussio n can be found in the Ap- 
pendix |A] Following iDavis et al.l ( 1201 ih . the starting point for the 
application of the GEV statistic is the cumulative distribution func- 
tion (CDF) 



PGEw(m) = Prob.()jj„,ax < m) 
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which gives the probability of finding a maximum halo mass m^^ 
smaller than m. In extreme value theory it has been shown that the 
CDF takes the following functional form ( Fisher & Tippett.192^) 



-PgevC") = exp-^- 
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where u = logigm is the random variable and a, f} and y are the 
shift, scale and shape parame ter of the distribution, respectively. 
As shown in lDavis et alj ( l20fll) . the parameters in the Poisson limit 
(on scales ^100 Mpc hr^ for the application to galaxy clusters) are 
found to be 



7 = n{> mo)K - 1, 
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where mo is the most likely maximum mass, n(> mo) is the co- 
moving number density of haloes more massive than mo, V is 
the comoving volume of interest and -t^I is the comoving mass 
function evaluated at mo. The most likely mass, mo, can be found 
toavis etalj20Ilh by perfoiming a root search on 
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Here = OmoPcrit is the mean matter density today, vo = 
[Sc/cr(mo, zjf' with primes denoting derivatives with respect to 
m a nd A, a and p are th e parameters of the Sheth-Tormen 
(ST; ISheth & Tomieiil ( Il999l) ) mass function, which we will use 
throughout the paper. 




Figure 1. Scheme for tiling the sky with subvolumes of size V, at a redshift 
z and subsequent measurement of the most massive cluster with mass in[^^^ 
in the volume. The faces of the subvolumes at redshift z are assumed to 
be squares with an angular extent S6 and all subvolumes have an extent in 
redshift space of 6z- 
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Figure 2. Number of square patches of angular scale SB (right axis of ordi- 
nate) as a function of angular scale for a survey area of 20 000 deg' (black 
line). In addition, the redshift dependence of a fixed comoving transverse 
distance between 100 — 300 Mpc ft"' (left axis of ordinate) is shown as a 
function of angular scale for the ACDM cosmology. For better orientation, 
the grey dotted lines denote the lines of 500 and 1000 patches, respectively. 



3 INTRODUCING THE CONCEPTUAL IDEA 

In order to observe the CDF, one has to sample a number of sub- 
volumes so that one may compile a sample big enough to recover 
the underlying distribution function. In the following, we study this 
approach in the framework of a hypothetical deep wide-field sur- 
vey covering an area of 20 000 deg^ and capable to detect clusters 
at redshifts above z = 1. 

The basic idea, as depicted in the scheme shown in Fig. [T] is to 
tile the sky with square patches of angular extent 66 (side length 
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of the patch) and depth Sz, which could be chosen to correspond to 
a given comoving length for a given cosmological model. In prin- 
ciple, these subvolumes could be placed on a spherical shell any- 
where in redshift space, but one has to keep several things in mind. 

(i) Naively we want to have as many patches on the sky as pos- 
sible, to get a better sampling of the underlying distribution. 

(ii) The subvolumes must have a minimum size, such that haloes 
contained in them can be considered to be uncorrelated and the 
Poisson limit is valid. Furthermore, larger volumes usually lead on 
average to a larger mmax, improving the detectability. 

(iii) The depth 6z has to be sufficiently large such that the red- 
shift determination of the most massive cluster is accurate enough 
to assign the object to the volume. 

(iv) The limiting survey mass should be low enough to allow to 
completely sample the peak of the distribution as shown in Fig. lAll 

(v) Moreover, it has to be ensured by the chosen redshift and 
selection function of the survey that one statistically finds a system 
in each subvolume. 

The number of patches as expected for a hypothetical survey area of 
20 000 deg^ is shown as a function of angular scale in Fig.|2] show- 
ing that in theory one expects 10^ - lO'' patches. But this number 
alone is meaningless unless one knows to what redshift a comov- 
ing transverse length-scale corresponds for a given angular scale. 
The refore , Fig. [2 ] also shows the comoving transverse distance (see 
e.g. lHogd j 



Dt = DnSe, 



(5) 



where Dc is the comoving (line-of-sight) distance, as a function 
of angular scale for five different length-scales between 100 - 
300 Mpc/;"'. One directly reads off that is possible to achieve at 
least a few hundred patches at rather low redshifts and up to several 
thousands for higher redshift and small length-scales. 
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Figure 3. Evolution of the equatioti-of-state parameter w{a) as a function 
of scale factor a. 
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4 GEV FOR DIFFERENT COSMOLOGICAL MODELS 
4.1 For a small fixed comoving volume 

The attempt of determining the underlying CDF distribution by 
measuring the mass of the most massive cluster in the subvolumes 
as outlined in the previous section leads naturally to the question 
what differences we would expect for such a distribution for differ- 
ent cosmological models. 

To address this question, we computed the GEV distributions 
for seven different cosmological models, comprising the fiducial 
ACDM model with (h, JIao, ^mo , os) = (0-7, 0-73, 027 , 0.81) based 
on the WMAP 1-year results jKomatsu et alj [201 i t), an inverse 
power law model (INVl), a super-gravity model (SUGRA), the 
normal and generalised Chaplygin gas, a topological defect model 
with w = -2/3 and an extreme phantom model with w = -3. 
A more detailed description of the used test cosmologies, includ- 
ing the respective linear overdensity thr esholds S^, as we ll as cor- 
responding references can be found in IPace et ah 1 teoiol) . For all 
models, we use the fiducial ACDM parameters given above with 
exception of erg which is scaled according to 

<5c.de(z = 0) 

0"8.DE = 7 TTT^'^s ■ 

Oc.ACDMtZ - ") 

jAbramo et al.ll2007l) . The evolution of the equation-of-state pa- 
rameter w(a) with scale factor a is shown in Fig. [3] for all mod- 
els except the phantom one. It should be noted at this point that 
the models mentioned above can be referred to as test cases (not 



Figure 4. Redshift evolution of the ratio with respect to the ACDM model 
for mo (upper-left panel), the GEV shape parameter, y, (upper-tight panel) 
and the scale parameter, /}, (lower-left panel) for six dilferent cosmologi- 
cal models, computed for a cubic volume with a comoving side length of 
lOOMpc/i-'. 



necessarily any more consistent with current observations), stand- 
ing for cosmologies with a substantially different geometric evolu- 
tion and/or history of structure growth. The fact that our method 
probes only the high-mass tail of the mass function allows to test 
for cosmologies that agree more or less on the background level 
with ACDM, but exhibit substantial differences for the (non-linear) 
growth history. 

As a starting point we compute the three GEV parameters a, /?, y 
and mo for a fixed cubic comoving volume V with a side length of 
L = IOOM pc/i~' placed at red shifts in the range of z e [0,2]. As 
discussed in lDavis et al.l ^201 ih . redshift evolution in n{> m) can be 
neglected for such a small scale, but it is still big enough to guaran- 
tee the validity of the Poisson limit (see also AppendixlAt. By using 
the same comoving volume for all models, we neglect for now the 
different evolution of the cosmic comoving volume V which en- 
ters in equation (|4) and therefore in all distribution parameters in 
equation In doing so, one lays emphasis on the influence of the 
different structure growth on the results. 

In Fig.|4]we show the results for y, /? and mg; a was left out since it 
basically coincides with ing. We show in Fig. |4]the relative devia- 
tion from the fiducial ACDM model as a function of redshift for the 
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Figure 5. PDFs for all seven cosmological models in four dilferent redshift 
intervals of extent as labelled, assuming an angular scale SO = 6 deg of each 
patch. 
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Figure 6. Redshift evolution of the ratio with respect to the ACDM model 
for mo (upper-left panel), the GEV shape parameter, y, (upper-right panel) 
and the scale parameter, /J, (lower-left panel) for six different cosmological 
models, computed for a patch with angular scale 68 = 6 deg and an extent 
in redshift space of Az = 0.5. The redshift z denotes the lower end of the 
volume. 
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different cosmological models. Since we normalized all the mod- 
els to today they all coincide at z = 0, but start to differ the more 
we go to higher redshifts. It is directly evident from the upper-left 
panel that mo is the most sensitive parameter leading to deviations 
with respect to ACDM of at least ~ 10% at z = 1 and up to ~ 80% 
(Chaplygin gas) at z = 2. The shape and scale parameters y and y6 
are much less affected and vary only by a few percent for the given 
redshift range which will in practise not be measurable. 
The first result is that one can hope to potentially distinguish differ- 
ent models (with a sufficiently different growth history) at higher 
redshifts via mo, or equivalently a. Of course, it will presumably be 
harder to control the uncertainties in the mass measurements when 
going to extremely high redshifts. 

4.2 Volumes with significant extent in redshift 

Towards an observationally more realistic case, compared to the 
small fixed cubic volume with L = lOOMpc/;"', three points have 
to be addressed first: 

(i) In order to assign clusters to the subvolumes, the redshift de- 
termination has to be precise enough. In the case of EUCLID this 



Figure 7. Evolution with angular scale 66 of mo (upper-left panel), the GEV 
shape parameter, y, (upper-light panel) and the scale parameter, /?, (lower- 
left panel) for seven different cosmological models in the 1.0 < z < 1.5 
interval. 



will be based on photometric redshifts. The accuracy is estimated to 
be given by cr- = 0.05( 1 + z) and keeping in mind that 100 Mpc A"' 
corresponds roughly to Az = 0.03 - 0.05 for the redshift range of 
interest, one must significantly extent the subvolumes in redshift 
space. 

(ii) An extension in redshift space, however, means that we can 
no longer neglect the redshift evolution of n(> m) within the vol- 
ume, since this would lead to a significant overestimation of w^ax 
in V. 

(iii) Considering that we want to study the effects of different 
cosmological models on the expected GEV distribution we have to 
take into account that for a fixed Az and patch size 69 the volume 
will be different for each model. Therefore, we expect a contribu- 
tion by the different evolution of the cosmic volume in the parame- 
ters of the distribution as well. 
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Figure 8. Sampled CDF (blue histogram) for the fiducial ACDM model (upper-row) and the INVl model (lower row), theoretical CDF from the GEV 
distribution and the fitted CDF (red line) for the redshift interval denoted in the upper-right and an angular patch size of i5f = 6 deg. We observed 500 patches 
and binned them in 75 sample bins. The numbers in the grey shaded area to the right give the theoretical and fitted values of the GEV ct-parameter. Only 
statistical errors are shown. 



Thus, it seems to be more practical to just fix the patch size 68 and 
a redshift interval Az instead of fixing the comoving side lengths of 
V. 

The remaining task is then to define an efi'ective comoving number 
density neff(> "i) in Az. For this, we compute the absolute number 
Nv(> m) of haloes more massive than m in the volume and define 
«eff(> m) = A'v(> in)/V. Of course it is no longer possible to use 
the expression from equation (|4), but one has to solve the following 
equation (see Appendix lAt 



dm 



dtrfi 



dm 



0, 



(7) 



where d iicf; / d m\,„^ is the efi'ective mass function evaluated at mo 
which is related to the effective number density ni:ff(> ™) via 



dm 



d;ieff(> m) 
dm 



(8) 



Now we can compute the GEV distribution for volumes with a sig- 
nificant extent in redshift space. In Fig. |5] the probability distribu- 
tion functions (PDFs) for all seven previously mentioned cosmo- 
logical models are shown for four different redshift intervals with 
Az = 0.5 placed at z = {0, 0.5, 1.0, 1.5) having an angular patch size 
of 56 = 6 deg. It can be seen that in the two low-redshift subvol- 
umes the difi'erences between the PDFs are rather small, apart from 
the very extreme phantom model which has a significantly differ- 
ent volume evolution. At high redshifts, however, the PDFs for all 
models start to significantly evolve away from each other. This dif- 
ference is caused by the difi'erent growth history and the difl"erent 
evolution of the cosmic volume for each of them. 
The redshift evolution of the ratio with respect to ACDM of mo, y 
and j3 is shown in Fig. |6] The important result is that the peak po- 
sition of the PDF shows the strongest difference with respect to the 
fiducial model (more than 10% for z > 1 and up to 50% for z > 1), 
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Figure 9. Impact of an uncertainty in the mass estimation on the sampled distribution in 1.0 < z < 1.5: from left to right we assumed a log-normal distributed 
mass uncertainty with o"„, =0.1, 0.2, 0.3 dex for the fiducial ACDM model (upper row) and the INVl model (lower row). 



whereas for y and /3 the differences are less than 10% and show only 
a weak redshift dependence. In Sect.|6]we will discuss another rea- 
son why 7 cannot be used to constrain cosmological models. This 
is a confirmation that trying to measure mo (or a) might be the most 
promising option. 

In Fig.|7]we show the same parameters but as a function of angular 
patch size 69 in the 1.0 < z < 1.5 interval, which also corresponds 
to a change of volume, leading to an increase of otq with 66, as 
expected. 



5 SAMPLING AND FITTING THE GEV DISTRIBUTION 

In order to address the question of how well one can really observe 
the cumulative distribution function of the most-massive haloes one 
has to create a simulated sample of the underlying CDF. An arbi- 
trary CDF can be sampled by means of the inverse sampling tech- 
nique, requiring an inversion of equation (|2) for which an analytical 
form exists. A value u = log[Q m is drawn by plugging in an equally 
distributed random number between and 1 into the inverted rela- 
tion. 

For the sampling procedure a survey-area of 20000 deg^ was as- 
sumed, divided into 500 square patches of 6 x 6 deg". We concen- 
trate on two high-redshift bins (1.0 < z < 1.5 and 1.5 < z < 2.0) in 
which we expect the biggest difference between the fiducial ACDM 
model and the test cosmologies. After the creation of a sample, 
which corresponds to the "observation" of the most massive ob- 
ject in each of the 500 patches, the CDF can be constructed in a 
straightforward way. First, one sorts the sample with the criteria of 



increasing mass, secondly one divides the mass interval into a num- 
ber of bins (75 in this case) and the final remaining task is to count 
the number of systems with a mass less-equal to the respective mass 
bin and divide by the total sample-size. We found the sample and 
bin numbers to be a good compromise of the stacking in the bins 
(see Sect. [6}, the sampling of the distribution and subvolumes big 
enough to find clusters in an observable mass range. The samples 
obtained in this way are shown for ACDM and the INVl model 
in Fig. [8] by the blue histogram where the black line denotes the 
theoretical CDF that has been sampled. 

Since in a real world application the black line is not known, 
we have to fit the observed CDF with the theoretical curve from 
equation ^ and determine the parameter values. From the results 
presented in the previous section it is obvious that only mg, or o- 
have the potential to be used for probing diflrerent cosmologies. The 
shape parameter, y, is far too noisy and the scale parameter, p, only 
weakly discriminates between the different models. Therefore, we 
fully concentrate on exploring what can be done by measuring a. 
The values o-tmc from the theory and the fitted value, ttfi,, for the 
ACDM and INVl models, are given on the right of each panel in 
Fig.[8]and similar plots in the following. The respective values are 
also depicted as circles on the black line for the theoretical CDF 
curve and on the red line for the fitted one, respectively. All GEV 
parameters of the sampled and fitted distributions shown in the fol- 
lowing Figs lSllll can be found in Tab.[l] 

The good news is that 500 patches in 75mass bins are suf- 
ficient to obtain the underlying value of a extremely well for all 
cases shown in Fig. [8] The statistical error of the fit is almost negli- 
gible and the systematic errors lo-truc - fffltl are small. The difference 
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Figure 10. Impact of an uncertainty in tlie mass estimation on the sampled distribution 1.5 ^ z < 2.0: from left to right we assumed a log-normal distributed 
mass unceilainty with cr„, =0.1, 0.2, 0.3 dex for the fiducial ACDM model (upper row) and the INVl model (lower row). 



between the two models shown in Fig. [8] is a multiple of the sys- 
tematic error, such that in this idealised case one would clearly be 
able to distinguish between the two models. The interesting ques- 
tion now is to study how well the procedure outlined above works 
as we commit significant errors in measuring the individual masses 
of the most massive objects, which is discussed in the following 
section. 



6 IMPACT OF UNCERTAINTY IN MASS ESTIMATES 

Unfortunately precise mass determination of galaxy clusters re- 
mains challenging, such that in order to understand whether the 
method proposed in this work is really applicable, it is necessary 
to study the impact of errors in the mass determination. To get an 
idea we model the error in the mass determination to be log-normal 
distributed with a cr„, = {0.1, 0.2, 0.3) dex. When drawing a value 
u - log[o m as discussed in Sect.|5]we change the value by adding 
a random error obeying the above mentioned log-normal distribu- 
tion. 

The results of this are shown in Fig. |9] for the redshift interval 
1.0 < z ^ 1.5 and the ACDM and the INVl model, where cr„, 
increases from left to right. We decided to show the results for the 
INVl model because it is among the models that show the biggest 
difference with respect to ACDM. As expected an increase in the 
mass uncertainty substantially alters the shape of the CDF, justi- 
fying the previous conclusion that the shape parameter, y, cannot 
be considered as a cosmological probe in this context. However, 
a, which is depicted by the red and black circles for the fitted and 
true value respectively, seems to be much less affected by the errors 



in the mass estimates. The fact that the sample is binned automati- 
cally implies a stacking of clusters with similar masses, which helps 
to reduce the scatter. Moreover, it should be mentioned that equa- 
tion l|2j delivers good fits in all cases even with substantial mass 
errors. 

As a cautionary advice it should be noted that, when following the 
procedure outlined above, one assumes the observed CDF to be the 
one of the true most massive clusters; however, observations will 
always provide the CDF of the most massive observed clusters in 
the patches. The difference between those two distributions stems 
from the fact that the most massive observed cluster might be a less 
massive cluster that scattered up due to errors in the mass determi- 
nation of the cluster. Incorporating this effect into the above anal- 
ysis is non-trivial and needs further study before the method could 
be applied to real data. To some degree the stacking might help to 
reduce the impact of this effect, but nevertheless the curves of the 
predicted observed CDF in Figs 1911 II should systematically shift 
to higher masses, since even if the mass of the true most massive 
cluster scatter s down it is inevitable that that a system with a sim- 
ilar mass scatters up. Thus the discerning power of the suggested 
method can only be used if the impact of mass eiTors is better un- 
derstood. 

A first inspection by eye shows that the difference o-fi, between the 
ACDM and the INVl model is substantial and gets less pronounced 
the bigger the mass uncertainty is. For the high-redshift sample 
1.5 < z < 2.0 depicted in Fig. [10] the situation is even better as 
the difference between the models gets more pronounced. 
As mentioned above it is crucial to have the limiting mass of the 
survey in the redshift interval of interest low enough to observe the 
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Figure 12. Difference between tfj^™^ and a™'''^' for all considered models. The grey shaded regions correspond to 1, 2 and 3 times the difference between 
the theoretical and the fitted value of a for the fiducial cosmology. From left to right the mass uncertainty increases from o"m = 0.1 dex to cr„, = 0.3 dex. The 
upper row shows the redshift interval of 1 .0 < z < 1.5 and the lower one the interval of 1 .5 ^ z < 2.0. 
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Figure 11. Impact of limiting survey mass of Mym - 10 Mq on the 
recovered CDF for the 1.0 < z < 1.5 interval for a lognormal distributed 
mass uncertainty with cr,„ = 0. 1 dex for the fiducial ACDM model (left- 
hand panel) and the INVl model (right-hand panel). 



low-mass tail of the CDF. From Fig. [TO]it can be inferred that a 
mass of roughly 10'"' Mq Ir' would serve even in the high-redshift 
case. Due to the robustness of the cumulative measure however, 
also a more conservative limiting mass of M\i„^ ^ lO'"*-^ A"' 
would still allow to recover the CDF almost equally well for the 
1.0 < z < 1.5 case as shown in Fig.[TT] Such a conservative value 
of Miim will be in reach of upcoming high-redshift surveys. The 
capability to discriminate between different cosmological models 
of the presented method is presented in Fig. [12] where the differ- 
ences between the value ffg, of the ACDM and the respective test 
cosmologies are shown. The panels display the differences with an 
increasing mass uncertainty from left to right, where the upper row 
presents the results for 1.0 ^ z < 1.5 case and the lower row for the 



1.5 < z < 2.0 case. The grey shaded areas show multiples of the 
systematic error lo-tme - Cfitl and we consider a model to be distin- 
guishable if its difference of cffl, is more than three times the sys- 
tematic error. It should be noted that this definition of detectability 
is extremely conservative, since the systematic error could be ac- 
counted for once the impact of the mass errors (lower mass clusters 
are accidentally mistaken as the most massive systems) is included. 
Once this is done only the statistical errors would remain which are 
much smaller and thus would lead to an substantial increase in the 
discerning power. 

As expected, the strongest ability to distinguish models is found for 
the lowest error in the mass estimates of cr,,, = 0. 1 dex and the high- 
redshift case as shown in the lower-left panel of Fig. [12] Here, all 
test models are significantly outside of the grey shaded area, which 
they enter in the panels for bigger cr„ to the right. In the case of 
cr,„ = 0.2 dex the models with the strongest difference with respect 
to ACDM could still be detected. 

From an observational point of view the 1.0 < z < 1.5 case how- 
ever is of bigger interest since the compilation of the required sam- 
ple will be easier. Also for this redshift range the o",„ = 0. 1 dex 
case is in principle able to distinguish all test models from ACDM. 
The more noisy cases could probably still be used as a first ACDM 
consistency check. Overall these findings show that the method dis- 
cussed in this work could be of great value for model-testing using 
a relatively small sample of objects, once the aforementioned treat- 
ment of mass errors is understood properly. 
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Table 1. Compilation of the GEV parameters of the sampled (true) and recovered (fit) distributions as shown in Figs lSIl l1 The values used in the aforementioned 
figures are arranged from top to bottom divided by horizontal lines where for the last two rows denoted with (-k), corresponding to Fig. ll II a limiting survey 
mass of Mii,^ = lO'"*-^ Mq instead of Mn^, = lO''' Mq ft"' was assumed. All fits are based on the initial values yo = -0.1, /?o = 0.1 and ag = 14.5. 
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-0.2852 ± 0.0074 


0.2196 ±0.0013 
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0.3 dex 
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14.4731 


-0.3067 ± 0.0078 


0.3123 ±0.0019 


14.4207 ± 0.0013 
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1.0 <z^ 1.5 


0.1 dex 
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0.0941 


14.5511 
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7 RESULTS AND SUMMARY 

In this work we propose a novel method for probing our cuiTent 
cosmological paradigm using the most massive galaxy clusters at 
high redshifts. The idea is to measure the cumulative distribution 
function (CDF) of the most massive haloes by tiling the sky with 
patches of the same area and identical depth in redshift space. We 
decided for a hypothetical survey of 20000 deg", divided into 500 
square patches of 6 x 6 deg" and Sz = 0.5 sampled in 75 bins, being 
a good compromise of the stacking in the bins, the sampling of the 
distribution and subvolumes big enough to find clusters in an ob- 
servable mass range. Such a survey could e.g. be the proposed EU- 
CLID mission or any other large-area survey able to detect galaxy 
clusters at high-redshifts with a sufficiently low limiting mass. All 
that is required is the ability to identify the most massive clusters 
in the individual subvolumes. 

The method that has been presented in this work has several ap- 
pealing advantages: 

(i) Using a cumulative measure is robust and the binning corre- 
sponds to a stacking of clusters of similar mass, reducing the impact 
of errors in the mass estimation. 

(ii) We are probing the high-mass tail of the mass function, 
well above the non-linear mass scale, where the difference between 
ACDM and alternative cosmological models is strongest. 

(iii) A limiting mass Mnm lO''' h^' is sufficient to sample 
the CDF even at high redshifts. For the 1.0 < z < 1.5 case also 
Mim - lO'"*-^ Mq /i"' would be sufficient to recover the CDF, which 
would definitely be achievable by future wide-area surveys. 

In order to measure the CDF, we fit a function obtained from 
generalised extreme value theory with three free parameters y, y6 
and a, which are the shape, scale and shift parameters, where the 
latter is closely related to the position of the peak of the underlying 
probability density function. Out of these three quantities the 
shift parameter, a, is the most sensitive one for model testing, 
whereas the shape and scale parameters are not sensitive enough 



and depend only weakly on redshift. Thus, we propose to use a as 
discerning measure between cosmological models. 
We show that by placing the patches at redshift z = 1.0 assuming 
an error in the mass estimation of cr„, = 0. 1 dex one can distinguish 
all models considered in this work clearly from the fiducial ACDM 
case. Such an accuracy of cr,,, = 0.1 dex in mass could be achieved 
by combining lensing^ and X-ray data. Going to higher redshifts, 
like z = 1.5 improves the discerning power further. Unfortunately, 
it is very doubtful that a sufficient mass accuracy can be achieved 
at such high redshifts by any currently available method. For larger 
errors in the mass estimation the discerning power is more and 
more reduced due the stronger impact of the convolution of the 
theoretical CDF with the distribution of the mass errors. This is 
also another reason why the shape parameter, y can not be used as 
discerning measure for different cosmological models. It should 
also be noted that by choosing the systematic error as measure for 
the discerning power we decided for a very conservative criterion 
which could, by correcting for the systematic error, be reduced to 
the much smaller statistical error. 

The main focus of this work was to present a novel method 
for probing the concordance cosmology using massive galaxy clus- 
ters. We demonstrated the potential of the method with a collection 
of toy-models, leading to encouraging results. The simplicity of 
the suggested scheme makes it in principle a perfect first test of 
ACDM, before a more demanding reconstruction of e.g. the halo 
mass function can be achieved. However, before this method can 
be applied in a real world scenario it is necessary to better under- 
stand the real impact of errors in the mass measurements on the 
observed sample. Particularly mistaking lower mass clusters as the 
true most massive system might cause a systematic shift to higher 
masses in the observed CDF undermining the potential as a cos- 
mological probe if not accounted for. Therefore this work can be 



For the 1 .5 ^ z < 2.0 case mass estimates from lensing can not be utilised. 
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considered as a first step to an applicable method for model testing. 
The next step will be to optimise the method further for application 
on cluster surveys and to better understand the impact of errors in 
the mass measurements as well as to study the discerning power in 
more detail for more realistic contenders of ACDM as well as to 
extend it to non-Gaussian models. 
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APPENDIX A: APPLYING GEV ON THE MOST MASSIVE 
HALOES 

In this appendix we outline the derivation of the most important re- 
lations for the application of GEV on the study of the most massiv e 
haloes in a volume of interest, as discussed in lDavis et al.l ( 1201 ll) . 
We start from the CDF given by 



(Al) 



/'gev('") = Prob.(mmax < ra) = PGEv(ramax) dm, 



In the following we will use both, the mass m and u = logj^ m as 
variables and decide case by case which is more convenient. The 
probability in equation dAlt to find a maximum halo mass smaller 
than m should equal the probability Po(m) to find no halo at all with 
a mass bigger than m. Thus the probably density function (PDF) 
PGEv(nj) is given by 



PoE\(m) = po(m) ■■ 



dPo 

dm 



(A2) 



If the volume of interest is large enough such that haloes can be 
considered to be unclustered, then Po(in) can be modelled by Pois- 
son statistic tDavis et al.,,201 li) 



Po(m) : 



A'' exp (-A) 
I] 



: exp [-«(> in)V] , 



(A3) 



where A = n(> m)V is the expected number of haloes more mas- 
sive than m, thus n{> m) is the comoving number density of haloes 
above mass m and V is the comoving volume. The number of oc- 
currences k is in the cuiTent case zero, leading to the last equality. 
The parameters cc,/? and y are obtained by performing a Taylor ex- 
pansion of the two CDFs and Pgev from the equations dBlt 
and l |A3t around the maxima of the distributions APqIAu and 
d/'cEv/dM, corresponding to the respective density functions. This 
gives 



Po(") = ^o(«o) + 
Pgew(u) = PaEwiiio) + 



dPoiu) 



d u 

d/'GEv(") 



d ;( 



{u-iio) + ... , 

(u - Uq) + . 



where the uo = logjo mo denotes the value of u at which the density 
functions have their maximum. Now, we equal the first two terms 
of the two expansions, the second order derivatives vanish by defi- 
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Figure Al. Probability density function of the GEV and Poisson statistics 
for a sphere of radius L = lOOMpc/i"' al z = as a function of u = 
log[Qm„,ax. The gray shaded area denotes the l-o" confidence region and 
the red line next to the peak (also known as mode of the distribution) is the 
expected value. 

nition and higher orders we neglect, hence one find;! 
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niQ In 10 

dm 



V = - 

•«a P 



1+7 



P 

("0 - a) 



(A4) 



(A5) 



where dn/dml,,,^ is the halo mass function evaluated at mo. To- 
gether with the relation for the mode (also known as the position of 
the peak) given by 



"0 ■ 



« + ^[(i+rr 

7 



1], 



(A6) 



one arrives finally at the equations for the three parameters 

(1 +y)(i+r) 



7 = n{> mo)V - 1, p - 



4^1 VmolnlO' 



a = logioOTo - -[(1 +7) ^ - 1] 



(A7) 



and p = 0.3. The analytic expression for nio is then found to be 

(All) 



mo 



a 1 ap(ayo )"<"""' 
2~27o~ 1 + (avo)-" 



+ ^ = 0. 



where vq is v evaluated at mo and primes denote derivatives with 
respect to mass. 



APPENDIX B: USEFUL RELATIONS 

In this section we summarise for the reader the most important re- 
lations for the GEV statistic as needed for this work. 



Cumulative distribution function (CDF) 



Pgem(u) = exp 



1+7 



-Uy 



(Bl) 



Probability density function (PDF) 



Pgev(") 



1 



1+7 



X exp 



1 + 7 



-i-i/r 



-i/r 



(B2) 



Mode - (most likely value) 
P. 



Mo = a + ^[(l+7)-''-l]. 

r 



Expected value 



Egev = a + -r(l - 7) . 

7 7 



Variance 



VARgev = ^ [f (1 - 27) - (1 - 7)] , 



(B3) 



(B4) 



(B5) 



where F denotes the Gamma function. All the equations given 
above are valid for 7 < 0, which is the case for the applications 
in this work. 



The most likely mass mo is obtained via the extremal condition 



d^Pn 



di<2 



, IdPo d^Po 

(In 10) m 1- m 

\ dm dm'^ 



which can be recast into 



dn d n I dn 
1- m- — - + mV\ - — 

dm d m'^ \ d m 



(A8) 



(A9) 



by using the definit ion of Pp from equation ( IA3t . When the mass 
function of Sheth & TormenI fl999) (ST) is incorporated in the 
above relation it i s possible to obtain an analytic relation for mo 
( iDavis et alj201 ij). The ST mass function is given by 



dn dlnv /ov _ 

-7- = — Ti A J— [1 +(ay) ''] e 

dm m- dlnm V 2;r 



(AlO) 



where p„, = flmoPciit is the mean matter density today, v = 
[6clcr{m,z)f' and the ST-parameters are A = 0.3222, a = 0.707 



^ When calculating dPo/dm one should keep in mind that dn(> m)/dm = 
—dn/dm. 



